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Abstract 

Climnlologiciil  vanalions  of  ihermal  and  halinc  buoyancy  tliixcs  are  investigated  in  the  Black  Sen.  Analyses  arc  pert'orined  to 
deieniiine  whether  or  not  thennal  buoyancy  flux  due  to  net  heat  flux  {or  haline  buoyancy  flux  due  to  freshwater  flux)  dominates  net 
buoyancy  tlux  in  the  Black  Sea.  Tlie  cfTect  of  the  two  types  of  buoyancy  tlux  are  examined  using  a  =^3.2  km  resolution  HYbrid 
C  oordinate  Ocean  Model  (IIYC'OM).  In  the  Black  Sea,  salinity  is  increased  by  Bosphorus  inflow  and  decreased  by  precipitation 
and  the  iiiHow  from  six  major  rivers.  Thus,  the  monthly  mean  discharge  values  from  six  major  rivers  are  used  as  additions  to  the 
precipitation  Held  in  the  model.  River  discharges  are  obtained  from  four  climatologies;  (1)  River  Discharge  (RivDIS),  (2)  Perry,  (3) 
University  Cor|xiriUion  for  Atmospheric  Research  (UCAR)  and  (4)  Naval  Research  Laboratory  (NRL).  Statistical  evaluations  of 
climatological  river  discharges  from  these  products  re.sult  in  similar  annual  mean  values.  However,  there  are  differences  in  the 
seasonal  cycle.  In  the  case  of  Danube,  which  has  the  largest  annual  river  discharge  contribution  of  ~6365  ni'  s  *,  RMS  differences 
for  river  discharge  \ allies  over  the  seasonal  cycle  are  within  ^2%  among  all  products.  The  Black  Sea  HYCOM  simulation  uses 
climatological  monthly  mean  atmospheric  forcing  (wind  and  thermal  forcing)  from  European  Centre  for  Medium-range  Weather 
Forecast  (ECMWF)  Re-analyses.  Buoyancy  Ilux  fields  obtained  from  the  HYCOM  simulation  demonstrate  that  thermal  buoyancy 
flux  dominates  haline  buoyancy  flux  in  all  months  except  March  when  the  basin-averaged  absolute  ratio  of  former  to  the  latter  is 
0.5  in  the  Black  Sea  On  the  contrary,  large  buoyancy  ratio  values  of  I  in  other  months  explain  the  buoyancy  is  much  more 
sensitive  to  variations  in  heating.  It  is  also  found  that  near  the  Bosphoais  Strait  in  the  Black  Sea  the  strongly  concentrated  source  of 
salty  w  ater  typically  penetrates  into  the  deeper  layers  as  a  plume.  Further  analyses  of  mean  buoyancy  fluxes  reveal  nonexistence  of 
deep  convection  in  the  Black  Sea  on  climatological  time  scales. 
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1.  Introduction 

The  ocean  is  dnven  by  a  combination  of  w  ind  stress 
from  the  atmosphere,  heat  fluxes,  freshwater  flux  resulting 
from  evaporation,  precipitation  and  ninotT  from  rivers. 
Among  tlicse,  the  freshwater  flux  is  of  particular  im¬ 
portance  since  the  upper  ocean  undergoes  a  regular  cycle 
of  convection  and  re-stratification  in  response  to  tlic  annual 
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cycle  of  buoyancy  lliixcs  at  the  sea  surface  (c.g.,  Webster. 
1994;  Marshall  and  Schott,  1999).  fherefore,  accurate 
dcicmiinaiion  of  freshwater  llu.xes  along  with  heat  lluxes 
at  the  sea  surface  is  of  particular  importance  for  many  types 
of  applications,  including  climate  studies  (c.g.,  Trenberth 
ct  al.,  2001). 

Freshwater  dischiirge  into  the  oeean  basins  has  become 
increasingly  important  in  global  climate  system  (Miller 
ct  al.,  1994;  W’ijnel.s,  2001;  Dai  and  Trenberth,  2002). 
fh'opcr  estimates  of  continental  freshwater  discharge  is 
particularly  essential  for  small  ocean  basins,  such  the 
Black  Sea  where  the  river  runofU' substantially  a  fleets  the 
salt  balance,  and  freshwater  lluxes  maintain  the  stable 
stratification  in  the  Black  Sea  (Ogu/  and  Besiktepe.  1 999; 
Stanev  ct  al  ,  2003).  In  addition,  there  are  many  un¬ 
certainties  in  the  existing  heat  flux  climatologies 
constnicicd  from  local  observations  which  arc  sparse  in 
both  time  and  space  (c.g.,  Sehnim  ct  al.,  2001). 

fhe  Black  Sea,  main  focus  of  this  paper,  can  be  con¬ 
sidered  as  an  estuarine  basin  as  there  arc  rivers  discharged 
in  to  the  continental  shelf  (Fig.  1).  In  the  northwc.stcm 
shelf  the  river  ninoff  is  maintained  by  Dimubc,  Dnieprand 
Dniestr.  hi  addition  to  these  rivers,  Rioni,  Sakar>'a,  Ki/i- 
lirmak  and  many  other  small  ones  discharge  into  the 
Black  Sea.  The  nanxnv  Bosphoms  Strait  plays  an  im¬ 
portant  role  in  detennining  the  ventilation  of  the  Black 
Sea  (O/soy  et  al..  2001 ).  Thus,  a  combination  of  heat  and 
freshwater  lluxes,  river  nmolTand  the  Bo.sphorus  outilow' 
can  have  substantial  impact  on  dynamical  features  of  the 
Black  Sea.  A  simple  examjile  is  that  the  amount  ofnver 
runofl'  in  the  continental  shelf  is  clo.sely  tied  to  the  Cold 


Intennediate  Layer  (CIL)  fonnation  (e.g.,  Ogu/  and 
Besiktepe,  1 999).  Buoyancy  due  to  river  ninotT  is  also  a 
contributing  factor  in  maintaining  the  basin-wide  cyclonic 
circulation  (Ogu/  et  al.,  1995).  This  becomes  very  im¬ 
portant  especially  in  the  northwcstcni  shelf  where  the 
majority  oflVcshwatcr  discharge  occurs.  In  particular,  the 
Danube  River  can  have  substantial  effect  on  the  ocean¬ 
ographic  features  of  the  western  Black  Sea.  This  riv  er  has 
monthly  mean  river  discharge  values  ranging  from  4(K)0 
to  9000  m^  s  '  on  climatological  time  scales  (\drosmarty 
ct  al.,  1997,  199S),  making  it  one  of  the  word's  largc.st 
rivers.  By  collecting  its  water  over  large  part  of  F'uropc, 
this  river  integrates  the  atmospheric  signal  and  creates  an 
amplified  forcing  in  the  western  Black  Sea. 

All  of  the  factors  mentioned  above  demonstnite  tlie 
impact  of  nver  mnolT  effects  in  the  Black  Sea.  Thus, 
accurate  modeling  of  salinity  is  of  particular  importance  in 
the  region.  Salinity  gradients  can  innuence  both  local 
circulation  dynamics  and  other  upper  occtin  fcauircs,  such 
as  the  vertical  stratification.  The  use  of  reliable  river 
discharge  values  is  therefore  essential  in  numerical  ocean 
modeling  studies  of  the  Black  Sea.  However,  the  problem 
is  to  decide  which  data  source  is  most  accurate  and  if  any 
arc  sufficiently  accurate.  Traditionally,  ocean  general  cir¬ 
culation  model  (OGCM)  .studies  in  the  Black  Sen  made  use 
of  the  river  discharge  values  from  two  .sources  (Altman  and 
Kumish,  19S6;  Staneva  and  Stanev,  199S).  Both  of  these 
data  sources  reported  river  mnolT  v  alues  constructed  from 
local  measurements.  Since  tlien,  there  have  been  other 
readily  available  climatological  river  discharge  values 
reported  from  a  few  global  data  sets.  One  aim  of  this  paper 
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1  ig  1 .  Miijor  nvers  discharged  into  tlic  Black  Sea.  Tlic  nioutli  of  each  river  is  sliown  with  numbers  from  1  to  6.  The  Ikispliorus  is  the  only  soiirvc  which  is 
considercil  as  negative  precipitation  (i  e..  evaporation)  source  in  the  mcxlel  simulation,  which  will  be  explained  in  Section  4,  to  close  the  evaporation 
precipitiition  budget  in  tJie  Black  Sea.  Kizihrmak.  also  known  as  Kizil.  is  the  longest  nv  er  in  Turkey,  ll  flow  s  «  1 150  km  from  the  central  l  urkish  plateau,  first 
to  the  southw  est,  then  to  the  north  and  nortlKast  into  the  iBack  Sea.  The  Sakatya  River  is  the  second  largest  river  discharged  into  the  Black  Sea  from  Tuikcy 
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is  to  discuss  monthly  and  annual  mean  ri\  cr  flow  values 
discharged  into  the  Black  Sea  as  constructed  from  these 
various  climatologies.  In  addition,  we  would  like  to  pro¬ 
vide  readily  a\ailable  river  discharge  values  which  can 
possibly  be  used  for  estuarine  and  OGCM  studies  in  the 
Black  Sea. 

Assuming  that  heat  and  freshwater  fluxes  from  opera¬ 
tional  models,  such  as  European  Centre  for  Medium-range 
Wcatlicr  Forecast  (ECMWF)  data  (Gibson  ct  al.,  1999), 
provide  complete  and  reliable  infonnation  at  all  temporal 
and  spatial  scales  in  comparison  to  climatologies  fonned 
from  local  observations  in  the  Black  Sea  (c.g.,  Schmm 
el  al.,  2001:  Kara  cl  al.,  2()05a),  they  can  be  used  for 
estimation  of  buoyancy  fluxes.  flowc\  er,  in  a  small  region 
like  the  Black  Sea,  the  grid  resolution  of  ECMWF 
(l.l25°x  1.125°)  is  not  appropriate  to  study  small  scale 
features  due  to  ri\  cr  runoff.  E\’cn  in  the  global  ocean,  there 
are  laigie  uncertainties  arising  from  spai^sc  obscrs^itional 
coverage  when  studying  the  river  ninoff  through  fresh¬ 
water  Iluxes  and  water  mass  lomiation  (Wijffcls  et  al., 
1992:  Doncy  ct  al.,  1998).  TTicrc  arc  also  problems  with 
air  sea  fluxes  from  model  outputs  (c.g.,  ECMWF), 
especially  near  the  coastal  regions.  This  is  due  mainly  to 
land  conuuTiination  of  atmospheric  variables  within  a  grid 
box  near  the  land  sea  boundary  (c.g.,  Kara  cl  al.,  2007). 
1  hcrcforc,  it  is  appropriate  to  calculate  sensible  and  latent 
heat  fluxes  based  on  the  model  SST  using  bulk 
foniiulations  rather  than  obtaining  them  directly  from 
ECMWF  (c.g.,  Kara  ct  al.,  2005b, c).  For  these  reasons,  an 
OGCM  can  be  useful  for  studying  the  impact  of  riverine 
and  buoyancy  iluxes,  especially  near  coastal  boundaries. 

fhcrc  are  various  studies,  focusing  on  mixing  processes 
at  vanous  regions  of  the  Black  Sea.  For  example,  Ogu/ 
cl  al.  (1990)  studied  internal  hydraulics  of  the  exchange 
processes  in  the  Bosphorus  Strait  from/to  the  Black  Sea 
based  on  a  simple  two-layer  numerical  model.  Ozsoy  ct  al. 
(2001 )  investigated  mixing  processes  near  the  Bosphorus 
Strait  and  some  coastal  boundaries  of  the  Black  Sea. 
Stane\  ct  al.  (2(K).^)  examined  the  Cold  Intennediate  Water 
(CIW)  1‘onnation  with  res[x:ci  to  buoyancy  fluxes  in  the 
iilack  Sea,  concluding  that  the  region  is  dominated  by 
dilution  of  suri'aec  waters  by  rivers  while  the  seasonal 
variability  is  controlled  by  the  air  sea  fluxes.  Koroiaex 
ct  al.  (2(X)6)  mentioned  about  possible  effects  of  the 
buoyancy  llux  on  the  Black  Sea  eireulation  dynamics. 
One  missing  part  in  these  earlier  studies  is  the  direct 
impact  of  ri\er  discharges  on  the  buoyancy  fluxes.  Es¬ 
sentially,  no  specific  consideration  was  given  to  river 
(lisehaQ^e  \  alues  especially  near  the  Bosphorus  Strait  and 
the  Danube  River. 

Fo  our  know  ledge,  fine  resolution  spatial  and  temponil 
variations  of  surface  buoyancy  fluxes  ha\  e  not  been 


examined  in  connection  wath  ri\  cr  discharges  in  the  Black 
Sea.  To  fill  in  this  gap,  the  main  focus  in  this  paper  is  to 
examine  freshwater  fluxes  with  a  particular  attention  to 
the  river  runoff  based  on  up  to  date  data  sources  and  then 
use  them  in  forcing  a  numerical  ocean  model.  We  will 
examine  heat  and  freshwater  fluxes  along  with  thermal 
and  halinc  buoyancy  fluxes  in  the  Black  Sea  using  a  fine 
resolution  (^3.2  km)  hybrid  coordinate  OGCM,  a  model 
which  includes  mixed  layer  physics,  ri\cr  forcing 
treatment  and  atmospheric  forcing. 

In  this  study,  we  will  analyze  climatological  mean  river 
discharge  values  obtainai  from  various  sources,  and 
investigate  differences  tunong  them.  Specifically,  we  put 
together  the  existing,  but  dispersed,  data  sets  and  assess 
them  carefully.  As  presented  in  this  paper,  in  detail,  these 
data  sets  are  also  made  a\ailablc  to  the  research  com¬ 
munity.  A  fine  resolution  OGCM  is  then  forced  witli 
optimal  monthly  river  flow^  values  in  addition  to  wind  and 
thermal  (c.g.,  solar  radiation,  near-surface  air  temperature) 
forcing,  ‘fhe  model  simulation  is  analyzed  to  study  spatial 
and  temporal  variations  of  buoyancy  fluxes  in  the  Black 
Sea.  Major  purposes  of  examining  the  buoyancy  fluxes  in 
the  Black  Sea  arc  to  (1)  detennine  the  relative  contribu¬ 
tions  of  heat  and  freshw  ater  fluxes  to  the  surface  buoyancy 
flux  in  the  Black  Sea,  (2)  identify  the  regions,  where  heat 
(freshwater)  flux  contribution  to  the  surface  buoyancy  llux 
is  significiUit. 

This  paper  is  oi^ani/ed  as  follows.  Section  2  intrtxiuccs 
the  flow^  values  for  the  rivers  dischaig;cd  into  the  Black 
Sea  from  various  climatologies  and  explains  how'  they 
were  constmeted.  Section  3  compares  monthly  mean  dis¬ 
charge  values  from  tlicsc  climatologies.  Section  4  de¬ 
scribes  OGCM  used  in  this  paper,  including  a  detailed 
examination  of  river  flow  treatment  in  the  model  and 
turbulence  pammctcrizations.  Section  5  presents  buoy¬ 
ancy  fluxes  obtained  from  the  model  simulation  over  the 
Black  Sea.  Section  6  examines  model-data  comparisons. 
Section  7  gi\  es  conclusions  of  the  paper. 

2.  River  discharge  climatologies  for  the  Black  Sea 

There  are  six  major  rivers  discharged  into  the  Black 
Sea  that  arc  examined  in  this  paper  (Table  1 ),  and  the 
location  of  each  river  is  show  n  in  Fig.  1.  Among  these 
rivers,  Danube  has  the  largest  upstream  area  of 
^^8()7,()(K)  kiiE,  followed  by  the  Dniepr  with  an  area  of 
^463,0(X)  km“.  Danube  is  the  major  European  river  to 
flow  from  west  to  cast.  It  originates  in  Germany  as  two 
smaller  rivers  called  Brigaeh  Breg.  These  join  in 
Donaucschingen,  flowing  south-east  for  a  distance  of 
about  2850  km,  to  the  Black  Sea  in  Romania  where  the 
Danube  Delta  is.  The  Dniepr  River  finds  its  source  in 
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Tabic  1 

A  list  of  rivers  discharged  into  the  Black  Sea 


River 

Station 

(‘ouniry 

(Lai,  Lon) 

Area  (kin‘^) 

Danube 

Power  plant 

Romania 

(45.2'  N,  29.7”L) 

S(),7000 

Dniepr 

Dniepr 

Ukraine 

(46.5°N, 

46  3()()0 

Rioni 

Sakochakidze 

Georgia 

(42.2°N.  4i.6^E) 

13,300 

Dnicslr 

Bendery 

I  krainc 

(4ri.2‘^N,  30.i°H) 

66,100 

Sakarya 

Botbasi 

Turkey 

(4l  FN. 

55,322 

Ki/ilirmak 

Ino/u 

Turkey 

(4L7°N,  36.0T4 

75,121 

The  country  represents  where  the  river  mouth  is  located,  i.atitude  is  the 
appriiximatc  latitude  of  the  river  mouth  and  longitude  is  the  approximate 
longitude  of  the  river  mouth,  both  of  which  arc  based  on  the  NRL  river 
climatology  as  desenbed  in  the  text.  Also  given  are  station  names  and 
their  upstrciim  areas  rcpoiled  from  the  RivDIS  climatology.  Note  that 
only  those  major  rivers  which  are  used  in  this  study  are  provided. 


noilh  Russia  and  tuns  south  eventually  discharging  into 
the  Black  Sea.  Dniepr  comes  from  Russia  through  Bclams 
and  then  Ukraine,  and  has  a  length  of  ~  2200  km.  Dniestr 
rises  in  the  Ukraine  and  flows  toward  the  Black  Sea.  Fora 
.short  distance  it  marks  the  border  of  Ukraine  and 
Moldovia. 

There  arc  several  monthly  or  annual  river  discharge 
climatologies  constructed  for  the  global  ocean.  These 
climatologies  also  include  river  flow  values  discharged 
into  the  Black  Sea,  one  of  the  focuses  of  this  paper.  Before 
examining  the  monthly  and  annual  mean  river  flow  dis¬ 
charge  values  constructed  from  these  sources,  a  brief 
explanation  about  each  data  source  is  prov'ided  here.  As 
described  below,  there  are  mainly  four  climatologies  (I 
througli  IV)  from  which  river  flow'  values  discharged  into 
the  Black  Sea  arc  obtained. 

(I)  River  Discharge  (RivDIS)  climatology  (Voros- 
marty  et  al.  1 997):  This  data  set  contains  monthly 
discharge  measurements  for  many  stations  located 
throughout  the  world,  fable  2  gives  el imato logical 
values  for  major  rivers  discharged  into  the  Black 
Sea.  Mean  monthly  discharge  in  m^  s  *  was  derived 
by  summing  all  available  discharge  measurements 
fora  particular  month  and  dividing  by  the  number  of 
measurements.  Annual  mean  values  were  then  cal¬ 
culated.  'fhe  montlily  mean  values  arc  shown  in 
f  ig.  2  for  the  rivers  discharged  into  the  Black  Sea. 
The  discharge  values  reported  in  this  data  set  arc 
measured  through  the  use  of  a  rating  curve  that 
relates  local  water  level  height  to  water  flow.  Tliis 
rating  curve  is  used  to  estimate  discharge  from  the 
observed  water  level.  The  rating  curves  arc  pe¬ 
riodically  rcchecked  and  re-calibrated  through  on¬ 
site  measurement  of  discharge  and  river  stage.  Site 
attributes  were  checked  for  consistency  through 


comparisons  with  the  United  Nations  Educational, 
Scientific,  and  Cultural  Organization  (UNESCO) 
published  scries.  Speeifically,  cheeks  were  made  on 
the  accuracy  of  the  site  names,  locations,  and 
contributing  drainage  areas. 

(ll)  Perry  climatology  (Perry  ct  al..  1 996):  Tlie  data  set 
provides  estimates  of  annual  mean  river  discharges 
for  98 1  of  the  largest  global  rivers.  It  is  meant  to 
describe  frcshw'ater  discharges  to  the  oceans. 
Discharge  values  for  as  many  rivers  as  possible 
were  gathered  from  as  many  sources  as  possible. 
Often,  this  meant  Uicre  were  several  values  for  a 
particular  river.  For  each  river,  obvious  errors  and 
duplicated  values  were  removed  as  explained  in 
PeiTy  et  a  I.  (1 996).  It  was  assumed  that  if  two  or 
more  sources  quote  the  same  v  alue,  then  they  arc 
probably  referencing  the  same  original  source,  for 
each  river,  the  mean  and  standard  deviation  of  the 
remaining  values  were  calculated.  Then  values  that 
were  more  than  two  standard  deviations  from  the 
mean  were  eliminated  as  outlying  values.  The  mean 
and  standard  deviation  were  tlien  calculated  a  sec¬ 
ond  time  and  those  arc  the  values  used  in  this  data 
set.  The  accuracy  of  the  measurements  is  not  giv  cn 
in  the  original  sources.  River  gauging  is  generally 
thought  to  have  an  accuracy  of  5  1 0%  but  the  actual 
accuracy  depends  significantly  on  local  conditions 
(Dingman  1 994). 

Tlie  annual  mean  flow'  values  for  the  Black  Sea  rivers 
were  compiled  from  different  sources  published  in 
the  literature  (Table  3).  In  this  list,  the  number  of 
sources  is  the  ones  used  to  dctcnninc  the  arithmetic 
mean  stream  How  for  a  given  river.  However,  this 
does  not  necessarily  equal  the  overall  number  of 
sources  found  for  that  river  since  duplicate  values 
and  outliers  arc  not  used  to  dctcnninc  the  arithmetic 
mean  stream  flow',  f  or  example,  tlicre  arc  1 4  different 
.sv^urccs  for  tlie  Danube  River  but  only  6  of  them  arc 
used  to  calculate  mean  How  value  (  fable  4). 

(Ill)  University  Corporation  for  Atmospheric  Research 
(UCAR)  climatology:  Anotlier  global  data  base 
from  UCAR  contains  real  time  monthly  averaged 
river  discharges.  Geographic  cov  erage  of  the  joint 
set  spans  all  major  rivers  and  oceanographic  basins 
around  the  global  ocean,  including  the  ones  dis¬ 
charged  into  the  Black  Sea.  This  data  set  contains 
monthly  river  discharge  rates  for  4425  locations 
around  the  world.  A  problem  w'ith  this  data  set  is 
that  tile  measurements  arc  often  too  far  from  tlie 
mouth  of  the  river  to  reliably  represent  the  discharge 
into  the  ocean,  limiting  its  uscftilncss  in  ocemi 
modeling  studies. 
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lablc  2 


Monthly  cliniiitological  river  di.scharge  values  (in^  s 

)  Irom  three  dilTcrent  climatological  data  sets:  (1)  RivDlS,  (2)  UCAR  and  (3)  NRI 

Month 

RivDlS  monthly  mean  climatological  river  discharge  values 

Danube 

Dniepr 

Rioni 

Dnicstr 

Sakaiya 

Kizilirmak 

Total 

Jan 

5940.7 

1369.0 

302.4 

207.1 

267.9 

212.9 

83(X).0 

Feb 

h2l9  5 

1598.7 

345.3 

293.7 

272.1 

255.4 

8984.7 

Mar 

7367  1 

1672.9 

429.8 

550.6 

295.2 

328.4 

10,M4.0 

Apr 

8574.0 

2477.6 

652.9 

615.1 

269.9 

308.2 

12,897.7 

May 

X937.9 

2893.1 

610.1 

460.1 

183.0 

231.2 

1.3,315.4 

Juii 

8315.7 

1616.6 

533,6 

502.8 

146.7 

157.0 

1  1,272.4 

.lul 

7122.5 

1057.6 

426.9 

475.5 

122.8 

118  2 

9323.5 

Aug 

5519.1 

941.9 

325.4 

.348.8 

1  10.8 

123.8 

7.369.8 

Sep 

4703.8 

841.5 

240.2 

288.4 

112.0 

147.4 

6333.3 

Oct 

4446.5 

979.8 

293,7 

247.3 

123.6 

167.6 

6258.5 

Nov 

4996.0 

1 1 1 1 .5 

356.6 

260.9 

218.4 

173.4 

7116.8 

Dec 

5S39.9 

1 240.9 

385.1 

250.2 

191.8 

202.6 

8110.5 

lA'AR,  month!}'  mean  climatological  river  i/ischargc 

'  values 

Jan 

5881.4 

1 369.0 

286.5 

187.7 

267.8 

212.9 

8205.3 

Feb 

6025.2 

1602.1 

350.9 

231.8 

272.3 

256.1 

8738.4 

Mar 

7255.5 

1672.9 

443.4 

486.5 

295.2 

.328.4 

10,481.9 

Apr 

8607. 9 

2477.6 

640.6 

587.7 

269.8 

308.3 

12,891.9 

May 

8898.9 

2893.1 

644.9 

401.5 

183.0 

231.3 

13,252.7 

Jun 

8185.5 

1616.6 

550.3 

382.9 

146.7 

157.0 

11,0.39.0 

.lul 

6917.8 

1057.6 

416.0 

386.8 

122.7 

117.8 

9018.7 

Aug 

5365.8 

941.9 

292.9 

310.1 

1 10.7 

123.8 

7145.2 

Sep 

459S.3 

841.5 

243.2 

245.3 

1 12.0 

147.4 

6187.7 

Oct 

4404.5 

979.8 

304.4 

219.4 

123.6 

167.6 

6199.3 

No\ 

4971.1 

1111.5 

321.3 

231.1 

218.3 

173.4 

7026.7 

Dee 

5850.9 

1240.9 

338.0 

220.9 

191.7 

202.6 

8045.0 

\RL  nuvtthly  meati  cltmafoloyical  river  discharge  values 

Jan 

5818.6 

1505.0 

303.3 

180.2 

301.9 

212.9 

8321.9 

Feb 

6091  6 

1757.6 

346.3 

255.5 

306.7 

255.4 

9013.1 

Mai 

7215.7 

1839.1 

431.1 

479.0 

.332.7 

328.4 

10,626.0 

Apr 

8397.7 

2723.8 

654.8 

535.1 

304.2 

308.2 

12,923.8 

May 

8754.2 

3 1 80.6 

61 1.9 

400.3 

206.3 

23 1 .2 

13,384.5 

Jun 

8144.7 

1777.2 

535.2 

437.4 

1 65.3 

157.0 

11,216.8 

Jul 

6976.1 

1162.7 

428.2 

413.7 

138.4 

118.2 

92.37.3 

Aug 

5405  6 

10.35.5 

.326.4 

30.3.5 

124.9 

123.8 

7319.7 

Sep 

4607. 1 

925.1 

240.9 

250.9 

126.2 

147.4 

6297.6 

Oct 

4355.1 

1077.2 

294.6 

215.2 

1 39.3 

167.6 

6249.0 

Nov 

4893.3 

1221.9 

357.6 

227.0 

246.2 

173.4 

7119.4 

Dec 

5719.8 

1364.2 

386.2 

217.7 

216.2 

202.6 

8106.7 

Total  rivei  tlischargcs  IVom  all  rivers  are  proMtlecl  by  month.  The  onginal  global  RivDlS  data  set  is  available  online  at  http: '  daac.oml.go\  rivdis  .  and 
the  global  IJCAR  data  set  is  available  at  http:  .  dss.iicar.edu  datasc{s/ds552.h  .  The  NRL  river  climatology  may  be  available  upon  request.  Monthly 
mean  discharge  values  for  Dniepr  arc  same  for  RivDlS  and  UCAR  except  Februaiy  as  both  data  sets  use  almost  identical  data  sources  during  the  same 
time  period  ( 1 952- 1 9X4). 


(IV)  Naval  Research  Laboratory  (NRL)  climatology 
(Barron  and  Sinedstacl,  2002):  The  hypotliesis 
used  in  constructing  the  global  NRL  river  data  set 
proposes  that  a  database  of  monthly  mean  river 
discharges  will  be  supeiior  to  a  database  of  annual 
means  in  its  ability  to  estimate  real  time  discharge. 
Tlie  NRL  river  data  set  eomes  from  Peny  ct  al. 
( 1 996)  which  had  one  mean  v  aluc  for  each  river  but 


the  set  was  converted  lo  monthly  values  to  be  used 
in  ocean  modeling  studies. 

A  summary  of  the  conversion  to  monthly  mean  values 
is  provided  here  for  the  rivers  discharged  into  the  Black 
Sea.  As  a  first  step,  monthly  stream  flow  numbers  were 
obtained  from  the  electronically  available  RivDlS  data  set 
(http:  Avww.daac.oml.gov/daacpages/fix'dis.html)  lo  get 
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I  ig,  2  C'iimalological  monihiy  inean  disciiarge  values  for  six  nvers  discliarged  into  ihe  Black  Sea  as  oblamed  from  Ihe  Rn  DlS  data  set 


most  global  (non-USA)  rivers.  Once  all  the  river  inflow 
values  were  obtained,  for  most  rivers  the  value  that  was 
largest  (near  the  river  outflow  )  and  the  most  eompletc 
(longest  time  scries  of  measurements  and/or  the  least 
amount  of  missing  values  and/or  the  most  recent  values) 
were  used.  Finally,  the  annual  mean  value  was  scaled  up 
the  value  reported  in  Peiry  ct  al.  (19%)  for  each  of  the 
RivDIS  rivers  if  the  latter  w'as  larger,  under  the  as¬ 
sumption  that  the  smaller  value  resulted  from  stations 
ilirther  upriver,  Asa  second  step,  monthly  mean  values 
were  processed.  The  pertinent  monthly  mean  data  values 
from  Riv  DIS  were  obtained  for  the  Black  Sea  riv  ers,  and 
each  monthly  average  relation  to  the  annual  mean  was 
figured  out.  All  of  these  relativ  e  v  alues  were  averaged  if 
there  was  more  than  one. 

For  river  stream  flow  or  discharge  data  sets  recording 
monthly  means,  which  average  to  an  annual  transport 
smaller  than  that  of  the  Perry  data,  there  arc  three  as¬ 
sumptions.  (1)  The  (larger)  Perry  mean  was  detennined 
using  data  downriver  of  the  station  used  to  derive  the 
monthly  means  and  is  therefore  a  better  estimate  of  the 

luble  3 


Number  of  sources  in  the  Perry-  d;il:i  sel  lo  be  useti  in  calculaling 
the  climatological  annual  mean  ri\er  flow  \  alues  discharged  into  the 
Black  Sea 


River 

Nunibcr 

of  souiees 

River  discharge 

Tola! 

Useil 

(nv  N  ') 

(km' s  ’) 

Danube 

14 

7 

6305.0 

203.7 

Dniepr 

10 

6 

1630.8 

52.2 

Rioni 

3 

3 

409.7 

13.1 

Dniestr 

4 

4 

326.3 

10.4 

S;ikai  Vii 

3 

3 

217.3 

7.0 

Ki/ilimi;ik 

1S0.5 

5.8 

Note  lliat  duplicated  sources  were  eliminalcd  from  total  sources  for 
annual  mean  discharge  calculaiions 


mean  annual  transport  at  the  river  mouth,  (2)  the  variability 
described  by  tlic  relative  size  of  each  river  among  the 
monthly  means  is  a  goixl  proxy  for  the  variability  at  the 
river  mouth,  and  (3)  thus,  tlic  montlily  means  at  the  river 
mouth  would  be  scaled  as  (dtita  set  monthly  mean)  ^  (Perry 
annual  mcan)/(data  set  annual  mean).  The  resulting  cli¬ 
matological  monthly  mean  v'alucs  are  show  n  in  1  able  2  for 
the  rivers  discharged  into  the  Black  Sea. 

3.  Black  Sea  river  climatology  comparisons 

The  river  discharge  climatology  is  fomied  for  each 
river  over  various  time  periods  based  on  the  av  ailability  of 
the  data  (  fable  5).  fhc  Peny  data  set  was  constructed 
using  annual  mean  at  individual  years  so  the  total  number 
of  years  represents  total  of  these  individual  years.  On  the 
contrary-,  we  constnicted  the  climatological  discharge 
values  from  RivDIS  and  UCAR  using  montlily  mean 
discharge  v  alues  for  each  consecutive  year  dui  ing  the  time 
period.  Annual  mean  discharge  values  shown  in  Fig,  3 
clearly  reveal  that  Danube  has  the  largest  riv  er  flow'  con¬ 
tribution  (=^70%). 

A  scatter  plot  of  river  di.scharge  v  alues  among  dilTcrcnt 
climatologies  clearly  demonstrates  good  agreement, 
while  there  exist  some  deviations  for  Dniestr  and  Sakary  a 
(Fig.  4).  Some  difTcrences  between  NRL  and  otlicr 
climatologies  arise  because  the  river  discharge  fioni  NRL 
is  scaled  to  the  value  proportional  to  that  from  the  Perry 
climatology.  Thus,  the  monthly  mean  value  is  a  linear 
operation  on  the  RivDIS  numbers  and  is  not  very  difi'erent 
from  them.  Specifically,  a  version  of  the  RivDIS  database 
w'as  used  as  one  of  several  inputs  to  the  NRL  database.  In 
the  case  of  the  rivers  discharged  into  the  Black  Sea,  most 
Riv  DIS  data  were  the  result  of  measurement  ending  as  far 
back  as  19K4.  Annual  mean  data  from  (Peny  el  al„  1996) 
w^as  more  recent  and  was  assumed  to  be  more  correct  and 
more  near  to  the  river  mouth.  Thus,  the  monthly  means  for 
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Inblc  4 

I  isi  ol  references  used  in  the  Perr^  data  to  obtain  annual  mean  Ilow 
values  lor  each  river  discharged  into  the  Black  Sea 


Discharge 

Reference 

Discharge 

Reference 

(m'  s“’) 

(iiF  s  ’) 

Danube 

Dniepr 

6550 

Milhnian  aiul 

1700 

Kanimerer  ( 19X9) 

Meade  (19S5) 

6500 

Kiimmcrei 

1700 

DyncMiis  and 

(I9S9) 

N  i  Is.son  ( 1 994 ) 

6500 

LM  S(  () 

1670 

van  der  Lceden  el  al. 

(l‘«5) 

(1990) 

6450 

Showers  (1979) 

1670 

Leopold  ( 1 962) 

64.50 

Dynesius  and 
Nil.sson  ( 1994) 

16(>() 

Kempc  el  al.  ( 199 1 ) 

6450 

(  >aya  ( 1 9S  1  ) 

1 660 

Degens  et  al.  (1 99 1 ) 

6450 

kempc  ct  al. 
tl991) 

1650 

Meybeck  (19SS) 

M.^0 

Mc>bcck 

tl9SS) 

1620 

Showers  ( 1979) 

6250 

Kempc  e(  al. 
(1991) 

1485 

LNI  SC'O  (19X5) 

62.50 

Degens  et  al 
(1991) 

1 160 

( Vaya  ( 19S 1 ) 

(>2.50 

van  del  1  eeden 
(1975) 

Rioni 

6200 

S/cstay  ( 19S2) 

409 

1  M.SCO  (19X5) 

6175 

\  an  tier  1  eeden 

420 

D>nesins  and 

ct  al.  (1990) 

Nilsson  ( 1 994) 

6175 

1  eopold  ( 1 962) 

400 

van  del  Lceden 
(1975) 

Sakar>'a 

Dniestr 

202 

van  del  1  eeden 
(1975) 

Ml 

LNi:SCO  (19X5) 

257 

Showers  ( 1 979) 

511 

Kempc  el  al.  (1 991 ) 

195 

1  \1  S( () 

(19X5) 

510 

Show  ers  (1979) 

Kizilinuak 

507 

van  dei  L eeden 
(1975) 

202 

1  \l.S(() 

(19X5) 

159 

v.in  del  1  eeden 
(L^s) 

Discharge  values  in  this  data  set  for  as  many  rivers  as  possible  were 
gathered  iVom  as  many  sources  as  possible. 


the  nvers  in  RivDlS  were  scaled  to  calculate  monthly  ^ 

means  proportional  to  the  annual  me<ui  by  multiplying  the 
Perry  annual  mean  and  then  dividing  by  the  RivDlS 
annual  mean. 

Statistical  analysis  is  pert'onned  to  further  examine  i 

errors  in  monthly  mean  discharge  values  between  the 
climatological  data  set  pairs.  Following  Muipliy  ( 
the  statistical  relationships  used  here  between  the  < 

monthly  mean  river  discharge  values  from  RivDlS  (A^ 
and  NRL  (>0  can  be  expressed  as  follows: 


ME  -  Y  -A', 


RMS 


n 


1  1,2 


t  I 


-  1 

NRMS^-  -  -  > 

,1 


/  I 


[X,  -  X){y,  y)/{axay). 

ft  *  ^ 


i  I 


ss  =  -  (^t/^.v)]"  -  [(r  -.v)/aAT  • 


(3) 

(4) 

(5) 


where  //“  12  (months),  ME  is  the  mean  error,  RMS  is  the 
root-mean -square  difTercncc,  NRMS  is  the  nonnalizcd 
RMS,  R  is  the  correlation  coefTicient,  SS  is  the  skill  score, 
and  A'(Y)  and  rrv  (rr>)  are  the  mean  and  standard 
deviations  of  the  RivDlS  (NRL)  discharge  values,  res¬ 
pectively.  Because  the  NRl^  climatology  is  based  on  an 
estimate,  it  is  always  considered  as  a  dependent  viu’iable 
(i.c..  F)  in  comparisons.  Note  that  RivDlS  is  considered 
as  a  reference  data  set  because  it  is  independent  of 
the  Perry  data  and  assumed  to  be  more  accurate  than 
UCAR  for  ocean  mc^eling  studies  as  the  IJCAR  data  set 
does  not  report  discharge  values  at  the  river  mouth. 
When  comparing  the  UCAR  climatology  with  the  NRL 


Tabic  5 


Time  pcnoils  dvcv  which  chmaiological  nvci  discharges  were  conslnicted 


River 

Perry 

RivDlS 

UCAR 

(“linratology 

Year 

Climatology 

Ycin 

('limatology 

Vear 

Danube 

1921  1994 

14 

1921  19S4 

64 

1921  2000 

80 

Dniepi 

1962  1994 

10 

1952-  1984 

33 

1952  1984 

35 

Rioni 

1975  1994 

5 

1965  1984 

20 

1928  1984 

57 

Dniestr 

1975  1991 

4 

1965  19X4 

20 

1881  1985 

76 

Sakarya 

1975  1991 

5 

1976-1985 

8 

1976  1985 

8 

Kizilirmak 

1975  1991 

2 

1976  1985 

8 

1976  1985 

8 

louil  number  of  years  for  the  chmatolog>  is  also  included.  I'hc  NRL  data  set  was  created  using  mainly  the  Perry  data  set  so  the  time  period  over  w  hieh 
the  NRL  climatology  was  constructed  is  the  same  as  the  Perry  data  set. 
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I  ig.  3.  Annual  nican  flow  values  from  various  climalological  dala  sources  for  Ihc  major  rivers  discharged  into  the  tflack  Sea.  I  he  annual  mean 
discharge  \'aliic  for  each  river  is  calculated  using  climatological  monthly  mean  discharge  values.  Also  provided  is  the  total  annual  mean  flow 
discharge  from  all  rivers.  Danube  and  Dniepr  provides  most  of  river  discharges  to  the  Black  Sea,  (>0.72”o  and  17.86”"o  of  the  total,  respectively. 
C’ontribiilions  to  the  total  river  discharge  from  other  rivers  arc  \‘cr)'  small  with  values  of  4.49"o.  3.57‘^'«,  2..3S'’  o  and  1 .9S”«  for  Rioni,  Dnicstr.  Sakarva 
and  Ki/ihimak.  respect i\ely.  The  percentage  values  were  calculated  using  annual  mean  values  from  the  Perry  climatology  as  shown  in  the  figure. 


climatology,  the  tbmicr  is  considered  as  the  reference  data 
set  (i.e.,  X)  because  it  is  ati  itidcpcndctitly  created  data  set. 

The  noti-ditnctisiotial  SS  which  iticlitdes  conditional 
atid  unconditional  biases  in  Eq.  (5)  is  also  used  because 
we  need  to  e.xatninc  more  than  the  shape  of  the  seasonal 
cycle  using  R.  The  SS  tneasures  the  accuracy"  of  NRL 
discharge  values  relative  to  RivDlS.  The  conditional  bias 
(^cond)  is  the  bias  in  standard  deviation  of  the  NRL 
discharges,  while  the  unconditional  bias  (/^uneond)  retleels 
the  mismatch  between  the  mean  NRL  and  RivDlS.  Skill 
score  is  1 .0  for  perfect  NRL  discharge  values.  1  he  NiLMS 
is  also  a  tneasiirc  ol  the  relative  distance  between  the 
e.stimate  and  reference.  The  scaling  factors  have  ditTcrent 
meanings  for  SS  and  NR  MS.  Both  non-dimensional 
metrics  are  used  to  obtain  further  infonnation  about 
differences  in  the  river  discharge  values. 

There  is  clearly  a  gex^d  agreement  between  Ri v  DIS  and 
IJC'AR  (sec  also  I  ig.  4),  with  small  ME  and  RMS  values 
for  each  river  (fable  6).  Almost  all  of  the  variance  (100%) 
in  the  RivDlS  di.scharge  values  is  reproduced  by  the  NRL 
and  UCAR  as  evident  the  square  of  all  R  values.  All  R 
values  are  >0.8.  For  the  1 2  monthly  river  discharge  values 
the  R  value  must  be  at  least  :r0.53  for  it  to  be  statistically 
difTerent  from  a  correlation  coefTieient  of  zero  based  on  the 
student's  t  test  at  the  95"  o  confidcnee  interval  (Ncter  et  af, 
1 98S).  I  he  SS  values  arc  very'  high,  being  very'  close  to  1 . 
This  indicates  almost  a  perfect  agreement  among  all  data 
sets  although  the  NRMS  values  for  UCAR  dischiirgc 

Nolo  here  llxit  iiccimicy  refers  to  the  nuteh  between  RivDlS  iiixl 
NRL  diseharge  vnliics,  nnd  skill  refers  to  the  NRL  disehiirge  accumey 
rclnlive  to  iho  Ri\  t)lS  chmalology 


values  vary  from  0.05"  o  to  10%  of  the  RivDlS  discharge 
values.  The  largest  enor  is  seen  for  Dnicstr.  In  the  ease  of 
RivDlS  N'crsus  NRL  and  UCAR  versus  NRL  pairs,  RMS 
and  Mli  values  are  slightly  larger  than  those  for  the 
RivDlS  versus  UC'AR  pair.  In  the  ca.se  of  Dnicstr  the 
NRMS  value  is  largest,  indicating  discrepancies  over  the 
annual  cycle  in  all  data  set.  A  remarkable  agreement  is 
ev  ident  at  Danube,  which  has  the  largest  discharge  values 
in  the  region.  The  river  discharge  estimate  for  Danube 
dilTcrs  only  with  ~2%  among  all  climatologies.  In 
summary,  as  evident  from  positive  and  large  SS  values, 
the  reasonable  agreement  in  discharge  values  for  all  riveiN 
is  quite  remarkable  given  the  uncertainties  and  measure¬ 
ment  errors  in  the  climatological  data  sets. 

It  is  emphasized  that  in  the  cases  of  riv  ers  discharged 
into  the  Black  Sea,  the  RivDlS  database  was  used  to 
detennine  NRL  river  database.  Most  Ri\  DIS  data  resulted 
from  measurements  ending  as  far  back  as  1984  as 
indicated  prev  iously  (  fahlc  5).  In  most  cases  the  annual 
mean  from  the  Perry  climatology  was  larger,  correspond¬ 
ing  to  a  closer  location  to  a  river  mouth.  The  exceptions  to 
this  w'crc  the  Danube  River,  v\'hcrc  the  mean  from  the 
Perry'  climatology  was  only  slightly  smaller,  ^2%,  and 
the  Dnicstr  River,  where  the  mean  was  ~  l3"o  smaller. 

4.  Ocean  model  description 

The  Black  Sea  simulation  is  perfonned  using  a  fine 
resolution  (~3.2  km)  HYbrid  Coordinate  Ocean  Model 
(HYCOM)  to  examine  buoyancy  fluxes.  Details  of 
HYCOM  equations  arc  given  in  Bleek  (2002).  The 
vertical  etx^rdinate  evaluation  for  HYCOM  is  discussed  in 
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•  RivDlS\ 

s  UCAR 

+  RivDIS  \ 

•s  NRI. 

0  UCARv: 

s  NRI. 

niscluirjic  (in  ^  s“’) 


l)iNcluirj:c  (in  ^  s”' ) 


1  iiz.  4.  Scatter  diagram  ot  climatological  monthly  mean  Mow  values  among  dittcrent  climatological  data  sets,  (‘omparisons  are  shtmn  at  each 
individual  river  discharged  into  the  tilack  Sea. 


C  liasstgiK'l  ct  al.  (2003);  the  vertical  mi.xing  algorithms 
arc  pro\  ided  in  ffalliwcll  (2004);  the  Black  Sea  model  set 
up  is  given  in  Kara  ct  al.  (20(k^b). 

The  model  combmes  the  advantages  of  the  iso- 
pycnal.  a  and  partial-cell  r-levcl  coordinates  within  a 
single  IVamework  in  optimally  simulating  coastal  and 
open-ocean  circulation  features  in  the  Black  Sea.  tn 
particular,  the  model  uses  hybrid  vertical  coordinate 
grid  generator  and  the  layered  continuity  equation  to 
make  a  dynamically  smooth  transition  from  isopycnal 
coordinates  in  the  stratified  ocean  to  a  terra  in-fol¬ 
lowing  coordinate  in  shallow  coastal  regions,  and  to  r- 
Icvcl  coordinates  in  the  mixed  layer  ancTor  unstratificd 
seas.  We  made  several  modifications  to  the  original 


hybrid  generator  routine.  In  summary,  the  original  ap¬ 
proach  can  lead  to  e.xccssivc  diffusion  when  remap¬ 
ping  layers  that  are  far  from  isopycnal.  Ihe  modified 
I  tYCOM  remapper  used  here  allows  the  profile  to  vary 
linearly  across  a  layer  when  the  layer  is  not  close  to 
being  isopycnal,  thus  significantly  reducing  numerical 
diffusion. 

4.  /.  Mixed  iaycr  para  meter  iza  lions 

Tile  assured  presence  of  unifomily  spaeetl  coordinate 
layers  in  the  upper  ocean  allows  the  fonnulation  of 
turbulent  near-surface  mixing.  The  K-Profilc-Paramc- 
tcrization  (KPP)  is  the  first  non-slab  mixed  layer  model 
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Stiilisliail  comparisons  of  Ihc  river  tlischaigc  valuer 

>  helwcon  vanous 

climalologieal  data  pairs 

for  each 

river  tlischargcd  inlo  Ihe  Black  Sea 

River 

RMS  (m"  s  *) 

Ml-;  (m’ s  ') 

R 

SS 

NMRS 

NMRS  ("'■>) 

Rix  DlS  ver.KtiS 

iCAR.  (lischaf};e  value.s 

Danuhe 

112  7 

S5.0 

1 .00 

0.99 

1.8C-2 

1.8 

Dniepr 

1.0 

0.3 

1.00 

1.00 

6.1c  4 

0.1 

Rioni 

2.V9 

5.S 

0.99 

0.96 

6.2c- 2 

6.2 

I3nicsli 

5S  2 

-50.7 

0.98 

0,81 

O.lc-0 

10.0 

Sakar>a 

0  1 

0.0 

1.00 

1.00 

4. Sc  4 

0.1 

Ki/ilinnak 

0.2 

00 

1 .00 

I  00 

1..3C-3 

0.1 

/^/r/)/.S  versa. ^ 

SRi  (liscltai'^e  values 

Danube 

137.0 

-  133.6 

1 .00 

0.99 

2.1e-2 

2.1 

Dniepr 

159,2 

147  4 

1,00 

0.93 

9.9C-2 

9,9 

Rioni 

1.3 

1.2 

1.00 

1.00 

3.0e  3 

0.3 

Dnicstr 

51.7 

48.8 

1.00 

0.85 

1.3c  1 

13.0 

Sakaiya 

260 

24.5 

1  00 

0,85 

1.3e-l 

13  0 

Ki/ilinnak 

0.0 

0.0 

1 .00 

1.00 

O.Oe  -  0 

0,00 

L  (\iR  ver.sus 

\  R/.  ilisehar^e  values 

Danube 

94.6 

48.6 

1.00 

0.99 

1.3c  2 

1.3 

Dniepr 

158  0 

147.1 

1.00 

0.94 

9.9c  2 

9.9 

Rioni 

24.1 

-7.0 

0.99 

0  96 

6,9e-2 

6  9 

Dnicslr 

24  6 

2.0 

0.98 

0.95 

6,3c  2 

6.3 

Sakaiya 

26.0 

24.5 

1.00 

0.88 

I..3C-1 

13.0 

Ki/ilinnak 

0  2 

0.0 

l.(M) 

1.00 

1.3c- 3 

0  1 

A  dclailed  description  of  ihc  siatisiical  parameters  is  given  in  the  icxl. 


(l  arge  ct  al..  1997)  included  in  HYCOM.  In  the  ocean 
interior,  the  contribution  of  background  internal  wave 
breaking,  shear  instability  mixing,  and  double  diiTusion 
(both  .salt  fingering  and  diffusive  instability)  arc  for¬ 
mulated.  In  the  surface  boundary  layer,  the  influences  of 
wind-driven  mixing,  surl'acc  buoyancy  fluxes,  and 
convective  instability  are  parameterized.  The  KPP  al¬ 
gorithm  also  accounts  for  the  influence  of  non-local 
mixing  of  temperature  and  salinity,  which  permits  the 
development  of  counter  gradient  fluxes. 

The  major  focus  of  this  paper  is  to  examine  surface 
buoyancy  lliixes  ov  er  the  Black  Sea.  Thus,  we  only  present 
parameterizations  which  arc  related  to  the  buoyancy 
lluxcs.  In  the  model,  buoyancy  flux  (B\)  and  buoyancy 
(li)  are  expressed  as  follows: 

B(  =  -  +  ^//(  r,.  .s; )( £■  -  P)s^,  (6) 

w 

mT^.s,)s.  (7) 


Po 


where  g  is  the  gravitational  acceleration  (9.K1  m  s  ^), 
«(7s,  Ss)  is  tile  tiicrmal  expansion  cocfTicicnt  *),  P 

S^)  is  salt  expansion  cocnicient  (psu  * ),  7’^,  is  the  sea 
surface  temperature  (®C),  5’^,  is  the  sea  surface  salinity 
(psu),  p  is  the  density  (kg  m  ^),/;,  is  the  reference  density 
(loot)  kg  m  ^),  is  the  net  heat  Ilux  at  the  .sc;i  surl'ace 
(\V  in  ^),  Cvv  is  the  ratio  of  interior  local  buoyancy  fre¬ 
quency  to  that  at  entrainment  depth  (a  constant  between  1 
and  2).  In  (6)  E  is  the  evaporation  due  to  latent  heat  (ms  * ), 
P  is  the  precipitation  due  to  rain  (or  snow  )  and  river 
discharges  (m  s  ’ ),  and  they  arc  discussed  in  Section  4.3, 
in  detail.  Expansion  coefficients  for  temperature  and  sa¬ 
linity  are  calculated  from  the  equation  of  .state  as  explained 
in  Brydon  et  al.  ( 1999). 

It  is  noted  that  in  (6)  includes  contributions  from  botii 
heat  flux  and  freslnvatcr  flux  at  the  sea  surface.  Positive 
(negative)  buoyancy  flux  indicates  a  buoyancy  loss  (gain). 
Surface  density  increases  (i.e.,  water  column  is  destabilized) 
if B^>  0,  and  surface  density  decreases  (i.e.,  w'ater  column  is 
stabilized)  ir^,  <0.  Net  heat  and  freshwater  fluxes  at  the  sea 
surface  arc  defined  as  positiv  c  quantities  into  the  ocean.  In 
addition,  B  in  (7)  is  a  function  of  the  relative  density  of  an 
ocean  parcel  compared  with  its  neighbor. 

4.2.  Black  Sea  moiiel 


^pjdS 

Po 


The  Black  Sea  model  has  a  1/25°  co.s(lat), 

( longitude X latitude)  square  mcrcator  grid.  Average  grid 
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resolution  is  3.2  kni,  ranging  from  a  minimum  of  3.05  km 
to  a  maximum  of  3.37  km  in  the  Black  Sea.  The  model 
includes  realistic  bottom  topography  constructed  from  the 
I  min  data  obtained  from  the  U.S.  Naval  Oceanographic 
Office  (NAVOCltANO).  The  climatologieally-foreed 
Black  Sea  IIYCOM  simulation  presented  in  this  paper 
was  performed  with  no  assimilation  of  any  oceanic  data 
including  SST.  This  is  done  to  examine  first  order  ctfeets 
of  atmospheric  forcing  on  the  model  simulation.  The 
model  is  initialized  using  the  temperature  and  salinity 
profiles  from  the  Modular  Ocean  Data  Assimilation  Sys¬ 
tem  (MODAS),  as  described  in  Kani  et  al.  (2()()5b).  The 
climatology  is  used  only  for  initialization  and  surface 
salinity  relaxation.  The  model  simulation  includes  no 
relaxation  to  any  data  except  a  relaxation  to  monthly  mean 
sea  surtacc  salinity  from  MODAS.  Tlic  model  simulation 
was  run  for  eight  years  using  the  6-hourly  climatological 
forcing  as  will  be  described  below.  It  takes  about  four 
model  years  for  a  simulation  to  reach  equilibrium.  For 
evaluation  of  the  model  results,  monthly  mean  model 
fields  were  fomied  from  daily  outputs  using  the  last  4 
model  years  (yciirs  5  through  8).  1'hc  reader  is  referred  to 
Kara  ct  al.  (2005b.c)  for  fiinher  details  of  the  model. 

BYCOM  reads  in  the  following  timc-vaiyring  atmo¬ 
spheric  forcing  fields:  wind  stress,  wind  speed  and  ther¬ 
mal  forcing  (air  temperature,  air  mixing  ratio,  shortwave 
radiation,  and  net  solar  radiation).  In  this  paper,  the 
IIYC'OM  simulation  u.scs  monthly  mean  wind/thennal 
forcing  based  on  E( "MWF  1 5  year  Re-Analysis  (HITA- 1 5) 
formed  over  1 979  1 993.  The  FRA  product  has  a  grid 
resolution  of  1 .  1 25°  x  1 .  1 25°.  Prc\  ions  simulations 
rcN'ealcd  that  the  use  of  the  liRA  forcing  results  in  reliable 
results  in  predicting  upper  ocean  variables  (Kara  ct  al., 
2005c).  The  model  run  is  perfonned  using  climatological 
monthly  mean  forcing  fields  inteqxdatcd  to  daily  values. 
A  high  frequency  component  (6  hourly)  is  added  to  the 
climatological  wind  forcing,  whose  details  and  applica¬ 
tions  are  given  in  Kara  et  al.  (2005 a). 

I  hc  net  heat  flux  absorbed  from  the  sea  surface  down 
to  depth  r  ((?(r))  is  given  as  follows: 

Q(Z]  -  C)(0)  f  (C^>i(0)  -  CToi(r)).  ( lO) 

0(0)  -CAw  4  Q\  (II) 

(?soi(r)/(.\>i(0)  (I  y)c\p(  r/0.5)  f  7cxp(  cA'par). 

(12) 

7  -  max(0.27,  0.695  5.7Ai»ak)  ( I3) 

where  ()(0)  is  net  heat  flux  absorbed  at  the  sea  surface, 
Qsoi(^)  ij'  shoi1\va\c  radiation  at  the  sea  surface. 


Q^,i(r)  is  remaining  (unabsorbed)  shortwave  radiation  at 
depth  r,  ()i  v,.  is  the  dow  nward  net  longwave  radiation, 
Q\  is  the  downw  ard  latent  heat  II ux  and  is  the  sensible 
heat  flux.  The  HYCOM’s  "surface"  heat  flux  is  not  (9(0), 
but  rather  the  near-surface  flux  absorbed  in  layer  I  {Q 
(3  m)  when  the  top  model  layer  is  3  m  thick).  Thus,  (^0) 
in  Eq.  (II)  does  not  include  None  of  ()soi(0)  is 

absorbed  at  the  surface  although  ( I  y)°  o  is  absorbed  vay 
near  the  surface.  The  model  reads  in  sjiatially  and 
temporally  varying  satellite-based  monthly  mean  attenua¬ 
tion  of  Photosynthetieally  Available  Radiation  (A>,\r) 
fields  to  include  effects  of  turbidity  in  the  shortwave 
radiation  penetration.  In  Eq.  (12),  the  red  penetration  scale 
is  2  m,  and  the  blue  penetration  scale  is  1  /Ap.\r. 

Net  solar  radiation  (net  shortwave  radiation  plus  net 
longwave  radiation)  is  so  dependent  on  cloudiness  that 
this  is  taken  directly  from  F:CMWF,  except  for  a  mod¬ 
ification  to  longwave  radiation  based  on  model  SST 
(Kara  ct  al.,  2()05c).  Both  Q\  and  in  F'q.  (II)  are 
calculated  using  model  SST  and  the  hulk  fonnulatic^ns 
that  use  stability-dependent  exchange  coefficients  at 
each  model  time  step  (Kara  ct  al.,  2()()5c).  leasing  Q\ 
and  (9s  on  the  model  SST  automatically  pro\  ides  a 
physically  realistic  tendency  towards  the  eorreet  model 
SST. 

4.3.  River  runoff  in  the  Black  Sea  nunlel 

H\'COM  treats  ri\crs  as  a  runoff  addition  to  the 
surface  precipitation  held.  The  flow  is  first  applied  to  a 
single  ocean  grid  point  and  smoothed  over  surrounding 
t'icean  grid  points,  yielding  a  contribution  to  precipita¬ 
tion  in  ms  ' .  In  the  Black  Sea  model,  there  are  a  total  of 
six  major  rivers  (Section  2). 

In  general,  the  Black  Sea  has  a  positive  freshwater 
flux  resulting  from  excess  precipitation  and  river  mnofl' 
over  evaporation.  This  positive  excess  can  be  con¬ 
sidered  to  be  equal  to  the  transport  in  the  Bosphoms 
Strait  (IJnluata  et  al..  1990).  The  reason  is  that  while  it 
is  not  thoroughly  confirmed,  the  Bosphorus  has 
generally  a  flow'  system,  such  that  less  saline  and 
lighter  water  exists  from  the  F3laek  Sea  via  Bosphorus 
at  the  surface,  while  the  saltier  and  denser  Mediterra¬ 
nean  w'atcr  flow  into  the  Black  Sea  as  undcrllow'.  fhus, 
the  transport  through  the  strait  eompensatc  for  the 
dilution  at  the  sea  surface  and  could  be  considered  to 
contribute  to  the  evaporation.  Because  the  simulations 
presented  in  this  paper  arc  carried  out  in  a  closed  basin 
eonUguration,  the  excess  of  precipitation  and  ri\cr 
l  unofTovcr  evaporation,  if  exists,  in  the  model  must  be 
removed  by  the  strait  outflow,  leading  to  a  zero  salt 
balance  in  the  Black  Sea. 
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Wc  express  the  net  water  balance  (/^nct)  in  the  Black 
Sea  as  follows: 

P  t  /^Kiver  t  ,  (14) 

where  E  is  eva])oration,  P  is  precipitation  due  to  rain  or 
snow;  i^  precipitation  due  to  ri\*crs,  and  is 

negative  precipitation  due  to  the  transport  from  the 
Bosphorus  Strait.  IlYC OM  uses  P  from  the  ERA- 1 5  as 
an  atmospheric  forcing,  and  /^Kncr  is  from  the  RivDIS 
which  is  comparable  to  NRL  and  UC"AR  climatologies 
(sec  Section  3). 

As  evident  from  the  net  w  ater  balance  gi\  en  in  Eq.  ( 1 4), 
the  Black  Sea  model  treats  rivers  as  a  “ninolT'  addition  to 
the  surface  precipiUition  field,  and  the  Bosphoms  Strait  as 
a  negati\  e  river  precipitation  (i.c.,  a  river  evaporation)  to 
clo.se  freshwater  flux  balance.  The  basin-averaged  values 
of  the  components  of  am  provided  in  Table  1,  Al- 
thougli  E  is  calculated  in  the  model  (i.c.,  it  is  not  a  pre- 
senbed  forcing)  using  simulated  latent  heat  flux  at  each 
time  step,  its  values  arc  provided  to  explain  freshwater 
balance.  Basin-averaged  annual  mean  E  (P)  values  are 
270.2  (221 .2)  km-  year  ‘  Ibr  FX'MWF. 

The  Pk,vvt  values  are  calculated  using  monthly  mean 
discharge  \  alues  from  Ri\  DIS  (sec  Table  2  for  individual 
monthly  discharge  values).  The  monthly  mean  Bosphorus 


Table  7 


liasin  averaged  iiionihly  and  aniuia!  mean  ovaporaiion  (/:),  precipitation 
iP)  v  alues  lor  Hie  Black  Sea 


HCMVVT 

Ri\or  How  tlischarges 

{ 6  Aiiuii) 

(knf  year 

(kni^  year  ') 

(knf  year  *) 

h 

P 

/V.vcr 

^  llosp. 

^  i«nal 

.Ian 

r  j 

1 

267.2 

260.3 

-215.5 

44.x 

46.5 

Keb 

-2A\.5 

220  3 

2X1.7 

-207.0 

74.7 

53.7 

Mar 

144.6 

1X5.3 

333.0 

190.1 

143.8 

184.5 

Apr 

-9]  1 

23X  4 

404  1 

’221. S 

1X2.3 

320.0 

May 

-ms.o 

231.4 

417.6 

-175.3 

242.3 

365.7 

.lun 

P)5.2 

233  0 

353.6 

1X5.0 

167.7 

205.5 

.lul 

-315.x 

222.6 

202.5 

-221.x 

70.7 

-22.5 

Aug 

-4064 

165.3 

231.2 

-234.5 

-3.3 

-244.4 

Sep 

4S0.S 

100.2 

10X.7 

240.3 

50.6 

344.2 

Oct 

426  0 

216  5 

1 06.2 

-245.1 

-4X.0 

-250.3 

\i)V 

-3X5.0 

.30X.0 

223.0 

-2.36  6 

-  13.6 

-01.5 

Dee 

281.4 

260.7 

254.2 

217.6 

Vv6 

24.0 

Alin 

2702 

221.2 

2X7.3 

216.7 

70.5 

0.5 

River  flow  diseliaiges  include  coiitnbutions  t'roni  the  total  of  ()  rivers 
(/Yivct)  i**^d  thmi  Bosplurnis  oulllow  ).NoIc  lhat  P^^  is  ihe  sum  of 

nver  discharges  useil  m  Ih  COM  (l  e..  ^  Ai<>sp )  The  values  in  knf 

year  '  i.s  obtained  by  nuiliipl>  ing  Ilow  values  in  km  ycitr  *  witli  the  surface 
area  of  the  Black  Sea  as  used  in  the  model  domain  (~4.()5  x  l(f  knT).  fhe 
unit  conversions  used  <ire  1  m  s  '  10  ^  .^0  ^  24  ^  60  nim  month  '. 

and  1  m  s  ’  ^  5(^5 ' 24  ■‘f>0^(i0  km  year  In  ihc  table  the  sign 

indicates  a  gam  to  ix'can,  and  sign  indicates  a  loss  from  txcan. 


outflow  values  are  taken  from  Staneva  and  Stancx  ( 1 90S). 
The  monthly  mean  Pr,vct>  and  Phos^  vtilucs  given  in  the 
table  arc  calculated  using  the  surface  area  of  the  Black  Sea 
(=s:4.05  10^  km")  to  be  consistent  with  the  basin-av- 
emged  E  and  P  v  alues.  The  annual  mean  bias  value  of 
0.5  km^  year  *  is  used  to  remove  the  excess  of  pre¬ 
cipitation  and  river  runoff  over  evaporation  in  the  model 
(T  able  7).  This  bias  mostly  accounts  for  runoff  \  alitcs  from 
many  small  rivers  wiiich  were  not  used  in  the  model 
simulations. 


5,  Buoyancy  fluxes  over  the  Black  Sea 


In  order  to  investigate  the  oceanic  response  to  the 
atmospheric  forcing  (heat  and  freshwater  lliixcs)  the 
buoyancy  fluxes  obtained  from  the  model  simulalion  arc 
presented  for  the  Black  Sea.  T  he  buoyancy  flux  expression 
given  in  Eq.  (7)  consists  of  two  icmis:  (I)  ihcmial 
buoyancy  flux,  and  (2)  haline  buoyancy  flux.  T  he  thermal 
buoyancy  flux  is  the  buoyancy  due  to  the  net  heat  flux  at 
the  sea  surface,  and  the  halinc  buoyancy  flux  is  the  one  due 
to  the  net  freshwater  flux  at  the  scii  surface,  fluis,  the  total 
buoyancy  is  rewritten  by  rearranging  (7)  as  follows: 


Bi 


Po 


+  ns.]. 


Thcnnal  hiK>ynncy  Ilux 


llahnc  biKnancv  tki\ 


(15) 


Hereinafter,  P  is  used  simplicity  to  represent  the  total 
precipitation  (i.c.,  the  sum  of  PR,ven  ^nd  P^osp  )•  4fic 
monthly  means  of  total  buoyancy  fluxes  along  with  its 
components  (thermal  and  halinc  buoyancy  fluxes)  arc 
presented  (big.  5).  The  monthly  fields  were  constnicled 
from  the  daily  model  outputs.  As  mentioned  earlier, 
negative  buoyancy  flux  indicates  a  buoyancy  gain  while 
positive  buoyancy  flux  indicates  a  buoyancy  loss.  The 
basin-averaged  monthly  mean  values  clearly  reveal  no 
strong  seasonal  cycle  in  buoyancy  fluxes  over  the  Black 
Sea  (big.  6).  A  buoyancy  gain  is  evident  in  most  of  the 
Black  Sea  from  March  through  September  on  the 
climatological  time  scales  (Eig.  5a),  Such  a  gain  .simply 
reveals  that  the  upper  ocean  is  effectively  sUibili/cd 
during  this  time  period.  This  buoyancy  gain  is  associated 
with  high  solar  heating  in  summer.  In  particular, 
relatively  large  (on  the  order  of  10  m"  s  ^) 
total  buoyancy  fluxes  arc  particularly  noted  from  May 
to  July.  On  the  other  hand,  there  is  a  substantial  buoyancy 
loss,  especially  in  the  northwestern  shelf  with  values 
>5x  10  ^  m^  s  '  in  January  and  even  10  '  in"  s  ^  in 
November. 
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(a)  I'olal  buoyancy 


JAN 


(b) 'riicrnial  buoyancy 


(c)  Halinc  buoyancy 


( 10  -«  in^  s  -^) 


I  ly.  5.  Spatial  and  temporal  variations  ot  buoyancy  fluxes  obtained  from  the  ==3.2  kni  resolution  IlYCOM  forced  with  hCMWl-  \\  ii)d  thennal  fluxes 
in  the  Black  Sea.  f  rom  top  to  bottom:  (a)  total  buoyancy  flux,  (b)  thennal  buoyancy  flux,  and  (c)  saline  buoyancy  flux  All  values  arc  in  10  ^  nr  s  \ 
and  the  results  are  shown  for  even-  other  month  The  surface  buoyancy  forcing  is  sum  of  the  thennal  buoyancy  flux  and  saline  buoyanc>  flux.  In  the 
color  bar  ^  denotes  smaller  (larger)  values  than  seen  at  the  edges 
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laii  Feb  Mar  Apr  May  Jim  Jiil  Aiig  Sep  Oel  Nov  Dee 
Fig.  6.  The  basin-averaged  monthly  mean  lota!  buoyancy  tliix  and  its  components,  They  are  based  on  Ihe  fields  given  in  (  ly 


Consistent  with  the  total  buoyancy  llux,  both  the 
ihcmial  and  ha  line  buoyancy  fluxes  tend  to  stabilize  the 
water  from  March  lliroiigh  September  over  the  most  of 
the  Black  Sea  (fig.  5b,c).  The  basin-averaged  buoyancy 
fluxes  further  con  firm  the  buoyancy  gain  in  these 
months  (fig.  6).  On  the  contrary,  both  thermal  and 
halinc  buoyancy  fluxes  result  in  destabilization  in 
January,  February,  November  and  December.  While 
annual  mean  net  heat  flux  is  zero  over  the  Black  Sea, 
there  is  no  reason  to  expect  basin -averaged  zero  mean 
annual  buoyancy  flux,  which  might  result  from  slight 
changes  in  mean  salinity  in  ver>'  long  time  scales. 

Using  Fq.  (15),  the  absolute  \  aliic  of  the  buoyancy  ratio 
([/^ )  of  the  tJiemial  and  halinc  buoyancy  flux  components  is 
given  by 


|A1  - 


y\T,.S.)Q, 


-  P}S, 


(16) 


where  the  ratio,  /?|,  indicates  the  relative  impact  on  upper 
ocean  buoyancy  of  heating  and  salinity  effects.  From  ( 1 6), 
;/?|  ^  1  (i.e.,  the  absolute  value  of  the  ratio  which  is  on  order 
of  unity)  explains  that  the  buoyancy  appears  to  be  equally 
alTeeted  by  heating  and  salinity  eflects.  In  other  words, 
heat  and  freshw^ater  fluxes  arc  of  the  same  magnitude.  In  a 
similar  analogy,  heat  flux  dominates  freshwater  flux  when 
R\ »  1  (i.e.,  the  buoyancy  is  due  mostly  to  net  heat  llux  at 
the  sea  surface),  and  freshwater  flux  dominates  heat  fluxes 
when  |/?|  1  (i.e.,  the  buoyancy  is  due  mostly  to  net 

freshwater  flux  at  the  sea  surface). 

The  spatially  and  temporally  \arying  \R\  values  in 
fig.  7a  calculated  using  thcmial  and  halinc  buoyancy 
flux\ allies  (fig.  5b,c)  demonstrale  that  net  heal  flux 
dominates  IVesh water  flux  on  the  Black  Sea  on  the 
climatological  time  scales.  Also  shown  arc  net  heat 
and  I'rcshwater  fluxes  at  the  sea  surface  (f  ig.  7b, e),  which 


are  used  for  calculating  thcnual  and  haline  buoyancy 
flux  values.  As  evident  from  the  figure,  there  are  large 
gradients  of  P-E  flux,  especially  near  Danube.  This  is 
because  river  discharge  values  in  HYCOM  are  considered 
as  additions  to  the  precipitation  fields.  The  same  is  also 
true  for  the  Bosphorus  which  is  treated  as  a  negative 
precipitation  source  in  the  model. 

Tlic  only  month  w^hen  thcmial  buoyancy  flux  destabi¬ 
lizes  the  upper  ocean  but  haline  buoyancy  flux  stabilizes  the 
upper  ocean  is  October  (Fig.  6).  Because  the  thermal 
buoyancy  flux  is  a  factor  of  16.8  larger  than  the  haline 
buoyancy  flux  in  October  (Fig.  8),  the  upper  cxrcan  de¬ 
stabilizes  and  turbulence  mixing  is  cxi^ectcd  during  the 
predawn  niinfall  in  tJic  Black  Sea.  The  net  surface  buoy¬ 
ancy  flux  also  destabilizes  the  winter  column  in  January, 
February,  Nov  ember  and  December,  which  is  an  indication 
of  the  generation  of  strong  turbulent  mixing.  Consequently, 
precipitation  tends  to  be  balanced  by  mixing  and  does  not 
produce  anomalous  surface  freshening.  The  freshwater 
fluxes  show  great  similarities  with  large  values  dunng  this 
time  [>eriod  (f  ig,  7b,c).  Thus,  turbulent  mixing  has  a  pro¬ 
nounced  monthly  variability  in  the  Black  Sea.  Such  a 
variability  in  the  turbulent  mixing  results  in  dissipation  of 
predawn  rainfall  that  mixed  downward  through  the  mixed 
layer.  In  HYCOM,  the  stability  of  stnitification  is  ensured 
by  the  fact  treating  Bosphorus  as  a  negative  precipitation 
source  is  very  local.  This  traitment  diKs  not  produce 
significant  mixing  over  large  areas  away  from  the 
Bosphoms  Strait  (see  Section  6). 

In  addition  to  the  thcnnal  buoyancy  flux,  an  exam¬ 
ination  of  the  freshwater  flux  clearly  reveals  that  preci¬ 
pitation  is  the  likely  reason  for  the  existence  of  shallow 
mixed  layer  formation  due  to  reduced  turbulent  mixing 
in  the  Black  Sea  during  spring  and  summer.  Note  that 
ver>'  shallow  mixed  layer  depths  were  noted  in  Kara 
et  al.  (20()5d).  For  example,  the  net  freshwater  flux  is 


A.li  Kara  cl  al.  Jiyurnal  of  Marine  Systems  '^4  (2iH)xS)  ‘^4  95 


SS 


JAN 


(a)  Buoyancy  flux  ratio 


(b)  Surface  heat  flux 


(c)  Surface  P  --  K  flux 
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Fig.  7.  (a)  Buoyanev  flux  ratio  (i.e..  the  ratio  ol'ihe  themial  buoyancy  tlux  to  saline  buoyancy  tliix),  (b)  net  heat  tliix  at  the  sea  surface  in  \V  m  ‘). 
and  (c)  net  freshwater  flux  at  the  .sea  surface  [K  /:  in  kg  m  ^  s  ’ ),  All  fields  arc  obtained  from  the  ^3.2  kni  resolution  1  lYC'OM.  Both  E  and  /^are  in 
m  s  and  they  arc  multiplied  by  the  density  of  ocean  water  at  the  sea  surtacc  to  obtain  the  flux  fonn  The  E  E  flux  values  in  the  figure  are  scaled  by 
to\  In  the  model,  we  adopted  the  con\ention  that  all  heat  flux  tenns  are  positive  into  the  ocean. 
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i  ly.  S.  I  he  basin-avera^cii  nionihly  mean  buoyancy  lalio  values  oblained  from  1!^  COM  in  ihe  Black  Sea.  Buoyancy  niiio  n  allies  aic  also  written 
above  hisiograms  for  each  monih.  Nolc  lhal  Ihe  buoyancy  ratio  (i.c.,  ihe  absolute  value  of  ihc  ratio  of  thcntial  buoyancy  flux  to  saline  buoyancy  Ilux) 
IS  iimtless. 


generally  positive  (i.e.,  into  the  ocean)  oxer  the  most 
of  the  Black  Sea,  and  this  is  particularly  seen  in  the 
western  Black  Sea  from  May  to  July  (1  ig.  7c).  In  these 
months,  having  P  £>0  implies  P>E.  Thus,  relatively 
large  prccipitatioti  values  dotiiinating  over  evaporation 
reduce  turbulent  mixing  and  produces  ihc  freshwater 
stratification.  In  a  prex  ious  study,  CTontn  and  McPhadeti 
(1999),  xvho  examined  the  relationship  between  rainfall 
atid  salinity  in  the  xvestem  equatorial  Paeiftc,  also  found 
that  rainfall  generally  produced  a  stable  shallow  layer 
when  the  net  surface  buoyancy  flux  stabilize  the  ocean 
suiTace  layer,  a  result  consistent  with  the  shallow  mixed 
layer  Ibniiation  in  the  Black  Sea.  In  particular,  the 
HYCOM  simulation  demonstrates  that  the  largest  and 
positix  e  net  freshxvater  flux  values  are  seen  from  March 
through  September  in  the  Black  Sea.  Such  a  feature  is 
consistent  with  estimations  of  net  freshwater  fluxes 
based  on  local  data  sets  (Sehmm  et  al.,  2001 ).  This  result 
further  confirms  the  existence  of  tlie  shallow  summer 
mixed  layer  due  partly  to  the  positive  freshwater  tlu.x. 
1  towever,  the  main  reason  for  the  shallow  mixed  layer  is 
the  net  heat  flux  since  the  thcmial  buoyancy  flux  is 
much  greater  than  haline  buoyancy  Ilux  during  summer 

(Tig  S). 

Finally,  a  possible  fonnation  of  deep  convection  in  the 
Black  Sea  is  investigated  by  examining  the  monthly  mean 
total  buoyancy  flux.  As  discussed  in  a  theoretical  and 
observational  review  study  by  Marshall  and  Schott 
(1999),  who  studied  formation  of  convection  oxer  the 
xarious  region  of  the  global  ocean,  whether  and  when  the 
deep  conxeetion  occurs  depends  on  the  seasonal  devel¬ 
opment  of  the  surface  buoyancy  Ilux  with  respect  to  the 
initial  stratification.  In  the  case  of  the  Black  Sea,  net  heat 
tlu.x  is,  for  example,  ^  100  W  m  ”  and  P-E  is  ~  1  m 
year  '  in  the  interior  during  summer.  This  implies  a 


buoyaney  tlux  of  ~  10  ^  nr  s  ^  ox  er  the  interior.  For 
stratification  typical  of  the  upper  regions  of  the  main 
thenniKlinc,  mixed  layers  do  not  reach  great  deptli  when 
exposed  to  buoyancy  loss  of  these  magnitudes,  perhaps  a 
several  hundred  meters  or  so.  Since  the  stratification  is 
very  strong,  and  the  total  buoyancy  tlux  is  not  sutficicnlly 
strong  (e.g.,  >  10  ^  m“  s~^)  in  the  interior,  convection  is 
not  expected  to  reach  greater  depths  in  the  Black  Sea  on 
monthly  time  scales.  In  other  xxords,  the  stratification 
inhibits  convection. 

The  existence  of  relatively  small  surface  thcmial  and 
haline  surface  fluxes  is  an  indication  of  the  fact  that  there  is 
almost  no  deep  conx  eetion  in  the  Black  Sea.  However,  the 
montlily  climatological  buoyancy  fluxes  may  expected  to 
be  much  xxxakcr  and  smoothed  than  the  episodic  buoy¬ 
ancy  forces  oj^rating  oxer  the  ocean  which  drive  the 
conxeetion.  The  conxeetion  is  an  instantaneous  process 
with  a  time  scale  of  few  days  but  I  lYC  OM  u.scs  monthly 
climatological  forcing  (i.c.,  longer  time  scale).  Thus,  there 
are  clctirly  dilTerences  in  time  scales  between  the  mixicl 
forcing  and  the  convection  episodes.  On  the  other  hand, 
xve  should  state  that  the  model  forcing  already  includes  6 
hourly  high  frequency  component.  Further  research  re¬ 
garding  possible  existence  of  convection  fonnation  on 
shorter  time  scales  (e.g.,  daily)  deserves  a  future  detailed 
study. 

6.  Model  validation 

This  section  presents  evaluation  of  model  results  and, 
in  particular  examines  whether  or  not  the  model  is  able  to 
reproduce  net  heat  Ilux  at  the  sea  surface,  which  was 
used  for  analyzing  the  buoyancy  flux  x'ariations  in  the 
Black  Sea.  Although  both  HYCOM  and  F'CMWF'  have 
zero  basin-averaged  annual  mean  heat  net  heat  Ilux  in  the 


A  B.  Kara  c(  al.  Journal  of  Xfarinc  Sysicm.s  "^4  (2()0S)  74  95 


^H) 


1  iji.  9.  The  biisin  iiveiagcd  monthl>  mean  climatological  net  tliix  al  the  sea  surface  from  FCMWF  over  1979  1993  and  ehmalologieally-foreed 
3.2  km  resolution  IIYC  OM  simulation. 


E3iack  Sea,  there  are  diOcrenees  especially  in  summer, 
Deeember  and  January  (Fig.  9).  The  basin-averaged 
RMS  net  Hu.x  diiTerenee  between  HYCOM  and 
FCMWF  is  24  W  in  "  o\er  the  seasonal  cycle.  This 
difference  e.xists  because  HYCOM  calculates  latent  and 
sensible  heal  tliixcs  based  on  the  model  SST  (rather  than 
FCMWF  SST)  at  each  time  step,  resulting  in  changes  in 
the  net  heat  II ux.  Such  differences  may  also  be  caused  by 
the  land-sea  mask  used  in  FCMWF  (Kara  et  al.,  2007). 
For  example,  atmospheric  forcing  fields  from  F:CMWF, 
as  used  in  the  1  lYCOM  simulation,  are  contaminated  by 
land  values  near  the  coast,  thereby  may  result  in 
unrealistic  SSTs. 

The  zero  net  heat  flux  in  HYCOM  is  maintained  using 
the  bulk  paramctcrizxitions  for  latent  and  sensible  heat 


fluxes.  Essentially,  tlicrc  is  a  feedback  between  SST  and  net 
fluxes.  Basing  fluxes  on  the  model  SST  automatically 
provides  a  physically  realistic  tendency  tow  ards  the  conect 
SST.  In  particular,  any  climatologically-forccd  (i.c., 
repeated  year)  simulation  will  have  zero  annual  net  heat 
flux,  w  ith  the  annual  mean  SST  automatically  adjusting  to 
make  this  happen.  Therefore  any  annual  mean  error  will  be 
in  SST,  rather  than  in  the  heat  flux  balance.  This  is  clearly 
evident  from  Fig.  10,  showing  the  basin-avcniged  monthly 
mean  SSTs  from  HYCOM  and  a  SST  climatology  as  well. 
The  monthly  mean-based  averaged  climatology  u.scd  for 
\*alidating  the  model  SST  is  based  on  the  NOAA^NASA 
Pathfinder  Advanced  Very  High  Resolution  Radiometer 
(AVHRR)  data  constnicted  over  1985  1997  (('ascy  and 
C  omillon,  1999).  This  data  set  represents  a  historical 


I  ig.  10  The  basin-averaged  climatological  monthly  mean  SST  from  the  Pathfinder  SST  climatology  and  the  3.2  km  resolution  climatologically- 
lorced  1 IVC  OM  simulalion  The  monthly  mean  Pathfinder  climatological  SST  used  in  this  paper  is  based  directl>  on  satellite  data  during  198.^  1997. 
It  has  a  resolution  of  -9  km.  I  his  climatology  is  based  directly  on  satellite  data  so  it  is  treated  as  the  “tmth'’  in  modekiata  comparisons.  The 
Pathfinder  SSf  climatology  was  interpolated  to  the  Black  Sea  HYC  OM  domain  before  the  basin-averaged  monthly  mean  SSTs  were  calculated. 
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(a)  Danube  zonal  section  at  45.2PN:  Climatological  mean  temp. 
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h)  Bosphorus  meridional  section  at  29.I4^K:  Climatological  mean  temp.  (®C) 
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Fig.  1 1.  ('ro.ss-scciii)iis  ot'climnlological  annual  mean  temperature  from  the  eddy-resolving  ~3.2  km  resolution  FIYC’OM  simulation:  (a)  Simulated 
zonal  seetions  in  I'Tont  of  the  Danube  Delta  al  45.2  l^N,  and  (b)  mendional  section  through  the  Bosphorus  Strait  at  29.I4°F..  The  thick  solid  line  in 
black  shows  mi.xcd  layer  depth  diagnosed  by  the  model  whose  caleiilation  is  explained  in  Kara  et  al,  (2(H)5d),  in  detail.  The  lines  inside  small  sub- 
panels  oflhc  both  ligures  illnstreites  the  locations  of  Danube  and  Bosphorus  seetions. 


reprocessing  of  the  entire  A\'1IRR  time  senes  using  con¬ 
sistent  SSI  algorithms,  improved  satellite  and  intcr-satcllitc 
calibration,  quality'  control  and  cloud  detection. 

Statistical  metrics  as  used  in  Section  3  are  also  used  lor 
assessing  the  accuracy  of  monthly  mean  HYCOM  SS  ls 
in  comparison  to  Pathfinder  SSTs.  The  model  gives  an 
annual  mean  bias  of -0.1 6  °C.  and  RMS  SST  difference 
of  0.98  °C  over  the  seasonal  cycle.  Almost  ^91%  of  the 
observ  ed  Pathfinder  SST  variance  can  be  c.xplaincd  by  the 


HYCOM  predictions  since  the  basin-averaged  correlation 
between  the  two  is  0.98. 

To  further  check  on  the  model  bias,  a  linear  regression 
analysis  is  performed.  If  it  is  assumed  that  the  HYCOM 
SST  would  fit  the  Pathfinder  SST  exactly,  then  a  fit  of  the 
observed  values  to  the  model  should  show  a  linear  slope 
with  a  zero  intercept.  However,  a  fit  of  the  observ  ed  data 
using  a  linear  fit  to  the  model  data  gives  a  relationship. 
Pathfinder  SST=  2.35  -  1 .  1 59  x  (H Y(X)M  SST)  ‘  error. 
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(;i)  Dunulic  /oiuil  section  at  45.21  N 
Climatological  mean  salinity  (psii) 


(b)  Bospliorii.s  meridional  .section  at  29.14 ‘K 
CMimatolo^ical  mean  .salinity  (psu) 
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I  'ig.  1 2.  The  same  as  I'm  1 1  hut  tor  cross-scction.s  of  climatological  annual  mean  .salinity  along  with  and  climatological  mean  salinities  in  March  and 
November 


The  linear  equation  can  be  reeasted  to  show  the  nature 
of  the  bias  as  follows.  Pathfinder  SST  =  1 4.78  t  1 . 1 59  x 
(HYCOM  SST- 14.78) terror,  where  the  annual  mean 
model  SST  is  14.78  The  eiTor  has  an  estimated  stan¬ 
dard  deviation  of  0.69  °C.  These  results  demonstrate  that 
HYCOM  is  biased,  with  the  bias  ehynging  at  about 
1 4.78  and  the  bias  increases  aw  ay  from  1 4.8  °C,  larger 
for  larger  values  and  smaller  for  smaller  values  of  the 
model  prediction.  If  1  lYCOM  has  only  random  sources  of 
error  in  prediction,  we  expect  the  slope  to  be  l.The  null 
hypothesis  is  tested  that  the  slope  in  the  model  is  one 
against  the  alternate  hypothesis  that  the  slope  is  not  1  .We 


found  a  /  value  of  0. 1 59/0.036~4.4  with  a  two  sided /> 
\  alue  ofO.OOl  6.  Thus,  the  null  hypothesis  is  rejected,  and 
it  is  accepted  that  HYCOM  has  non-random  sources  of 
eiTor. 

HYCOM  is  also  validated  using  subsurface  tempera¬ 
ture  and  salinity  fields.  Major  features  of  thelhcnnohaline 
fields  in  the  Black  Sea  am  the  diluted  surt'ace  water 
(salinity,  which  is  about  two  times  smaller  than  in  the 
Mediterranean  Sea)  and  the  CIL,  a  layer  where  the  tem¬ 
perature  shows  a  persistent  minimum  of  <8°C.  Tlic  latter 
generally  extends  to  depths  50  100  m,  where  the 
downward  pmpagation  of  the  themial  signals  from  sea 
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siirlacc  IS  blocked  b>  ihc  cxircmely  stable  strati ficalion. 

I  his  vertical  stratification  viirics  locally  as  a  consequence 
ol  dynamics  and  specific  distribution  of  sources. 

H  YC'OM  presented  in  this  paper  realistically  simulates 
most  of  the  above  mentioned  features,  for  example, 

1  ig.  1 1  shows  vertical  cross-sections  of  climatological 
annual  mean  of  subsurface  temperatures  through  two 
importiint  regions  where  the  model  Black  Sea  is  diluted 
(the  area  of  Danube  River  plume)  or  receives  salt  {the  area 
of  the  Bosphorus  Strait).  These  Helds  arc  calculated 
by  averaging  the  individual  daily  model  output  during 
4  years  of  integration.  The  w  arm  surface  layer  seen  in  the 
annual  mean  temperature  e.xtcnds  to  only  about  20  40  m. 
Below'  this  depth  the  cold  water  persists  throughout  the 
year.  The  local  minimum  gives  an  indication  about  CIW 
having  its  source  on  the  continental  slope.  This  cold 
plume  is  trapped  by  the  slopped  bottom  and  pmpagates 
southwards  along  the  western  coast,  clearly  revealing  the 
IT^'COM  success  in  simulating  coastal  processes. 

The  annual  mean  salinity  pattern  illustrates  the  fact 
that  the  major  dilution  of  the  surface  layers  is  caused  by 
the  river  runoff  (1  ig  12).  The  compensation  of  the 
excess  of  freshwater  mer  most  of  the  F^fack  sea  surface 
is  prescribed  in  the  model  in  a  narrow  \  icinity  of  the 
Fiosphorus  Strait.  The  strongly  concentrated  source  of 
positive  buoyancy  penetrates  into  the  deeper  layers  as  a 
vertical  buoyant  plume.  The  plume  is  vertical  down  to 
about  50  m,  w  hich  is  approximately  the  depth  of  the 
shelf  north  of  the  Bosphonis  Strait.  Below'  this  depth  the 
cuiTcnt  behaves  as  a  gra\  ity  flow'  reaching  depths  more 
than  100  m  before  losing  its  identity  (fig.  12b),  which  is 
also  mentioned  by  Korolaev  et  al.  (2001 ).  The  transport 
by  the  buoyant  plume  pro\  ides  a  continuous  source  of 
salty  water  helping  to  maintain  the  stratification  of  the 
Black  Sea.  The  phenomenology  of  the  simulated  out¬ 
flow  justifies  the  approach  chosen  in  this  paper  to  in¬ 
troduce  the  strait  outflow  as  a  negative  precipitation  at 
the  sea  surface. 

I  hc  sensitivity  of  salinity  in  the  upper  layer  to  the 
temporal  variability  in  the  freshwater  is  also  revealed  by 
an  examination  of  subsuiTace  salinity  fields  (1  ig.  12a,b). 
The  major  differences  are  t>bscrved  close  to  the  source 
areas  (the  Danube  Della  and  the  Bosphonis  Strait).  Fferc, 
one  important  result  has  to  be  made  clear.  Although  the 
chosen  periods  (March  and  November)  concspoiul  to 
periods  of  relatively  higher  ( lower)  precipitation  (ninoff), 
the  salinity  in  the  suifacc  layer  is  larger  in  March.  This  is 
explained  by  the  larger  mixing  in  the  winter  period  and 
can  be  seen  in  the  more  difTusiv  c  halocline  in  March.  On 
the  opposite,  the  seasonal  thennoclmc  reduces  mixing  in 
the  upper  layer  making  the  signature  of  frcslnvater  more 
clear  in  fall.  The  volume  occupied  by  the  fresh  coastal 


current  is  highly  variable  throughout  the  year  demonstrat¬ 
ing  that  this  buoyancy  signal  shapes  to  a  large  extent  the 
dynamics  of  the  Western  Black  Sea. 

7.  Summary  and  conclusions 

An  examination  of  freshwater  fluxes  along  w'iih  net 
heat  fluxes  at  the  sea  surface  through  the  buoyancy  flux  is 
essential  in  the  Black  Sea  because  they  influence  sim¬ 
ulation  of  upper  ocean  quantities.  Similarly,  an  accurate 
modeling  of  buoyancy  fluxes  is  of  particular  imjiorlance 
in  the  region  due  to  relatively  large  sal  inity  gradients  in  the 
vicinity  of  riverine  or  estuarine  dischmge.  Given  that  the 
river  ninolT  has  large  influence  on  the  dynamics  of  the 
continental  shelf  by  changing  the  heat  and  freshwater 
fluxes,  especially  in  the  northwestern  shelf,  a  database  of 
accunite  river  flow  estimates  is  useftil  to  support  v  arious 
types  of  studies,  including  observ'ational  and  modeling 
cfTorls  in  the  Black  Sea.  For  e.xamplc,  use  of  the  reliable 
river  database  also  enables  an  OGCM  to  more  realistically 
represent  v  ariability  and  di.stribution  of  near-shore  salinity 
in  the  vicinity  of  large  rivers.  It  also  prov  ides  a  centralized 
source  for  estimates  of  river  discharge  of  all  riv  ers,  ev  en 
smaller  rivers  which  might  be  insignificant  in  a  coarse 
resolution  OGCM  but  might  be  significant  in  a  high- 
resolution  OGCM.  in  this  paper,  a  specific  attention  is 
therefore  giv  en  to  rivers  discharged  into  the  Black  Sea. 

River  runofl  values  discharged  into  the  Black  Sea  am 
obtained  and  compared  to  each  other  from  four  climato¬ 
logical  data  sources.  Based  on  tlie  statistical  analysis, 
monthly  mean  river  discharge  values  from  all  data  sets  arc 
usually  close  to  each  other,  and  this  is  true  for  all  six  major 
rivers  (Danube,  Dniepr,  Dnicstr,  Rioni,  Sakaiya  and 
Kizilimiak).  Altliougli  the  montlily  mciui  river  discharge 
values  for  Danube  and  Dnicstr  show  some  differences 
when  comparing  RivDIS,  UCAR  and  NRL  climatologies, 
tlicse  differences  aie  negligible.  In  tlie  ca.sc  of  the  Danube 
River,  w'hich  has  the  largest  discharge  in  comparison  to 
other  rivers,  the  monthly  mean  discharge  values  from 
RivDIS  are  only  1.8%  (2.l®o)  different  from  the  UCAR 
(NRL)  climatologies.  Similarly,  the  discharge  values 
from  UCAR  arc  ^  f  .3%  different  from  RivDfS,  indicating 
a  remarkable  agreement  between  the  two  data  sets 
even  tlioiigh  UCAR  reports  discharge  values  far  away 
from  the  river  mouth.  Using  these  river  discharge  value  we 
also  analyze  freshwater  flux  balance  in  the  Black  Sea. 
Overall,  the  rivei's  provide  ^  287  km^  year  ^  precipi union 
accounts  for  ~220  knr^  year  \  evaporation  removes 
^270  knv^  year  and  outflow  through  the  Bosphoms  is 
^217  km^  year  ^ 

Following  the  analysis  of  the  river  data  sources  and 
freshwater  llux  balance,  a  climatologically-forced  Black 
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Sea  simulation  was  pcifonncd  to  calculate  net  buoyancy 
flux  along  with  its  components,  i.e.,  buoyancy  due  to  net 
heat  llux  (thennal  buoyancy  11  ux)  and  buoyancy  due  to 
freshwater  tlux  at  the  sea  surface  (ha line  buoyancy  flux). 
The  ftne  resolution  (^3.2  km)  HYCOM  simulation  is 
aimed  to  discuss  the  relative  contributions  of  heat  and 
salt  Iluxcs  at  the  ocean  surface  to  the  surface  buoyancy 
flux  m  the  Black  Sea.  A  buoyancy  gain  (i.c.,  the  upper 
ocean  stabilizes)  is  ev  ident  in  the  Black  Sea  from  March 
through  September.  On  the  contraiy,  a  buoyancy  loss 
(i.c.,  the  upper  ocean  destabilizes)  in  other  months.  The 
thermal  buoyancy  flux  is  much  greater  than  the  halinc 
buoyancy  llux  in  all  months  except  lor  March  when  the 
freshwater  flux  dominates  the  net  heat  tlux.  Thus,  it  is 
concluded  that  the  seasonal  variations  of  the  heat  flux  is 
the  major  contributor  to  the  buoyancy  llux  in  the  Ellack 
Sea.  October  is  found  to  be  the  only  time  when  thermal 
and  halinc  buoyancy  Iluxcs  have  different  signs.  In 
particular,  thermal  buoyancy  flux  destabilizes  the  upper 
ocean,  while  halinc  buoyancy  llux  stabilizes  the  upper 
ocean  in  October.  Because  the  thermal  buoyancy  flux  is 
much  larger  than  the  halinc  buoyancy  llux,  turbulence 
mixing  is  expected  during  the  predawn  rainfall.  While 
the  thermal  flux  is  seasonally  more  dominant,  local 
cITects  of  halinc  contribution  can  be  very  important  at 
some  regions.  We  also  speculate  that  since  the 
magnitude  of  total  buoyancy  flux  is  not  sufficiently 
strong  in  the  Black  Sea,  convection  should  not  be 
expected  to  reach  greater  depths  in  the  Black  Sea. 

HYCOM  is  able  to  predict  net  surface  heat  Iluxcs 
with  a  RMS  difference  of  24  W  m  ‘  in  comparison  to 
those  from  HC'MWF  over  the  seasonal  cycle.  IIYC  ()M 
also  maintains  net  zero  net  heat  flux  at  the  sea  surface  by 
using  bulk  parametenzations  to  calculate  latent  and 
sensible  heat  Iluxcs  based  on  the  model  SSTat  each  time 
step.  C’omparisons  of  basin-averaged  monthly  mean 
IIYCOM  SSTs  with  those  from  a  satellite-based 
Pathfinder  data  set  demonstrates  that  the  model  gives 
an  annual  mean  SS  f  bias  of-0.l6  °0,  and  RMS  SS  f 
dilTcrcnce  of  0.98  '^C'  o\  cr  the  seasonal  cycle.  Thus,  the 
Black  Sea  HYCOM  is  able  to  predict  both  net  heat 
Iluxcs  and  SST  with  reasonable  accuracies. 

fhe  use  of  river  forcing  as  an  addition  to  the  pre¬ 
cipitation  field  in  ocean  models  is  fairly  often  used, 
although  sometimes  it  is  handled  as  a  spread-out  land- 
amoff  c\enly  distributed  at  the  coastline  rather  than 
rivers.  There  is  not  much  justification  for  it  locally  (i.c., 
for  representing  where  river  water  advects  to),  but  in 
principle,  our  experience  demonstrates  that  it  improves 
the  basin-wide  salt  balance.  In  addition,  while  the  non¬ 
existence  of  typical  deep  convection  is  revealed  in  this 
paper  on  climatological  time  scales,  a  detailed  study  is 


needed  to  further  examine  possible  foniiation  of  sueh 
events  on  various  time  scales. 
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